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We study the crossover of the quantum Loschmidt echo (or fidelity) from the golden rule regime to 
the perturbation-independent exponential decay regime by using the kicked top model. It is shown 
that the deviation of the perturbation-independent decay of the averaged fidelity from the Lyapunov 
decay results from quantum fluctuations in individual fidelity, which are caused by the coherence 
in the initial coherent states. With an averaging procedure suppressing the quantum fluctuations 
effectively, the perturbation-independent decay is found to be close to the Lyapunov decay. We also 
show that the Fourier transform of the fidelity is determined directly by the initial state and the 
eigenstates of the Floquet operators of the two classically chaotic systems concerned. The absolute 
value part and the phase part of the Fourier transform of the fidelity are found to be divided into 
several correlated parts, which is a manifestation of the coherence of the initial coherent state. In 
the whole crossover region, some important properties of the fidelity, such as the exponent of its 
exponential decay and the short initial time within which the fidelity almost does not change, are 
found to be closely related to properties of the central part of its Fourier transform. 

PACS number(s): 05.45.Mt, 05.45. Pq, 42.50.Md, 76.60.Lz 



I. INTRODUCTION 

The quantum Loschmidt echo, M(t), measures the 
overlap of the evolution of the same initial state under 
two slightly different Hamiltonians in the classical limit, 

M(t) = |($ |ex P ( lj fft)exp(-^ t)|*o}| 2 - (1) 

Here Hq is the Hamiltonian of a classically chaotic system 
and 

H = H + kV, (2) 

with k being a small quantity. 

This quantity characterizes the stability of quantum 
motion under a small change of the Hamiltonian, e.g. by 
the interaction with the environment, it is thus called "fi- 
delity" and of great interest in the fast developing field 
of quantum information [1,2]. 

Moreover, a relationship between this quantity and the 
Lyapunov exponent that characterizes the classical chaos 
has been established analytically by using the semiclassi- 
cal theory by Jalabert and Pastawski [3] , which has been 
confirmed numerically in several models [4-9]. There- 
fore, this quantity has attracted a great attention from 
the community of quantum chaos [10-13], 

It has been shown that several regimes exist. In the 
golden rule regime of the perturbation parameter, above 
a perturbative border, the fidelity has been found to have 
a simple decaying behavior, M(t) oc exp(— Y^t) [4,6,14], 
where Yt, is the half width of the local density of states 
(LDOS) pl, which has a Lorentzian form (Brcit-Wigner 
form), 



with Yl = 2-KpU 2 . Here p is the density of states and U 
is the typical transition matrix element of kV between 
eigenstates of Hq. As the perturbation parameter k is 
increased, the exponential decay of M(t) deviates from 
exp(— Yr,t), when Yt, becomes comparable to the band 
width of Hq and pr, deviates from the Lorentzian form 
notably. A decay with the Lyapunov exponent A of the 
underlying classical chaotic dynamics, M(t) oc exp(— At), 
can appear in this regime of the perturbation parame- 
ter, when the (effective) Planck constant is small enough 
and the initial state is chosen suitably, e.g., a narrow 
wavepacket [3]. An interesting feature of this decay of 
M(t) is that it is perturbation-independent, in the sense 
that it is irrelevant to the strength of the perturbation 
kV and is determined by the classical behavior of the sys- 
tem Hq (and also that of H due to the smallness of the 
perturbation in the classical limit). Based on the semi- 
classical theory and the random matrix theory analysis, 
a transition of M(t) from the Pi-decay to the A-decay 
has been conjectured to occur at Yt, = X [4]. Numerical 
results in several models support the conjecture [4,7,9]. 



The validity of results of the scmiclassical theory de- 
pends on the value of the (effective) Planck constant. 
When the (effective) Planck constant is not small enough, 
the validity should be checked by direct quantum me- 
chanical calculations. In fact, the behavior of the fi- 
delity in this case is still not quite clear. For example, 
whether there is a sharp transition from the TL-decay to 
the A-decay, or the transition occur in a finite regime of 
the perturbation parameter. It is even unclear whether 
M(t) could decay with the Lyapunov exponent. Indeed, a 
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perturbation-independent but slower than exp(— At) de- 
cay has been observed in the kicked top model [4], the 
mechanism of which is still not clear. Meanwhile, the ran- 
dom matrix theory seems unsuitable for such problems, 
since perturbation-independent decay of M(t) is usually 
initial-state-dependent, the feature of which is hard to be 
captured by the random matrix theory treatment. 

The quantity used in this paper in analyzing properties 
of the fidelity M(t) is its Fourier transform, denoted by 
Fm{E) in what follows. As will be shown in Sect. II that 
F M {E) is determined directly by properties of the ini- 
tial state and of the (quasi)energy eigenstates of the two 
classically chaotic systems concerned. The role played 
by Fm{E) in understanding the behavior of M(t) is sim- 
ilar to that by the LDOS for the survival probability 
P(t) = \{a\cxp(-iHt)\a)\ 2 . 

In this paper, properties of Fm{E) and their relations 
to those of M(t) will be investigated numerically in the 
kicked top model [15]. The fidelity in this model has been 
studied in Rcfs. [4,14,13] and here we will concentrate on 
issues unaddressed in the previous work. Specifically, in 
addition to properties of M(t) related to the function 
Fm(E), we focus ourself on the following problems: 

(i) The mechanism of the perturbation-independent 
decay slower than exp(— Xt). 

(ii) Due to the small difference between H and H in 
the classical limit, the classical counterpart of the fidelity 
M(t) should change quite slowly when t is small enough 
(see Ref. [9] for numerical results in another model). 
In the quantum mechanical similar phenomenon 
should be of practical interest in quantum computing. 
However, it was not clear whether or not the phenomenon 
exists. 

The paper is organized as the follows. The function 
Fm is introduced in Sect. II. In Sect. Ill, we shall present 
our numerical investigation of the problems mentioned 
above, in particular, for properties of the amplitude and 
the phase of the function Fm(E) with initial coherent 
states. The perturbation-independent decay slower than 
exp(— Xt) will be shown to be due to large quantum fluc- 
tuations in individual fidelity, which are caused by the 
coherence possessed the initial coherent states, when the 
(effective) Planck constant is not small enough. It will 
be shown that the fidelity M(t) does not decay obviously 
within a short initial time, the length of which is deter- 
mined by some properties of the function Fm(E). The 
relationship between the decaying rate of M(t) after the 
short initial time and some properties of the function 
Fm{E) will also be shown numerically. Conclusions and 
discussions will be given in Sect. IV. 

II. FUNCTION F M (E) AS THE FOURIER 
TRANSFORM OF M(T) 

In this paper, the eigenstates of the two Hamiltonians 
Ho and H are denoted by \a) and respectively, with 
eigenenergies E a and Ep, 



H \a)=E a \a), H\/3) = E p \0). (4) 

The expanding coefficient of \(3) in \a) is indicated by 
Cf) a , C fja = (a\/3). The LDOS p L (E) of an eigenstate 
| a) of Hq is determined directly by the eigensolutions 
the two Hamiltonians concerned, 

p L (E) = J2\CfJ a \ 2 S(E-(E f ,-E a )). (5) 



As is known, when the initial state |$o) is an eigenstate 
of H 0} the fidelity amplitude m(t), 

m(t) = ($ |cxp(«#t)cxpHff t)|$o), (6) 

is just the Fourier transform of the LDOS pl(E), and 
the form of the LDOS is important in the study of M(t), 
which is just the survival probability P(t) in this case. 
However, as shown in Ref. [8], the LDOS can not explain 
the perturbation- independent decay of M(t). In fact, the 
perturbation-independent decay of M(t) appears only for 
some special set of initial states, that is, it is initial-state- 
dependent. Although a rigorous condition is still lacking, 
a sufficient condition for such initial states is expressed 
in a rough way that an initial state should be a narrow 
wavepacket [3]. When the decay of M(t) is initial-state- 
dependent, one should use a quantity that is more general 
than the LDOS in analyzing properties of M(t). 

The function Fm, the Fourier transform of the fidelity 
M(t), is a suitable candidate, which is introduced as the 
follows. Let us first express the fidelity amplitude m(t) 
in the form 

m(t) = J f m (e)e- iet de, (7) 

where 

,f m (e)=J2^S[e-(E a -E p )}, (8) 

a,/3 

A a& = <$ |/?>M*o>qSa- (9) 

Note that the quantity A a p defined in Eq. (9) is invariant 
under the transformation 

|a>->e*«|a>, |/?}^e^|/3>, |$ ) - e**"> |$„>, (10) 

where Lp a ,ip/3, and ipo are arbitrary phases. That is, it 
is irrelevant to the relative phases among the states |a), 
and |$o}- The fidelity M{t) can be expressed as 

M(t) = \m(t)\ 2 = J F M (E)e- lEt dE, (11) 

where 

F M (E)= J f m {e)t m {e-E)de. (12) 
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The quantity Fm (E) is in fact the correlation function of 
/ m (e). Note that Fm{E) can be calculated directly from 
properties of the states |a), |/3), and the initial state |$o}- 
In the following section, we will investigate numerically 
the relationship between properties of M(t) and those of 
Fm{E). To this end, it is convenient to write Fm{E) as 

F M (E) = p F (E)e ie ^ E l (13) 

Some arguments can be given to the reason why the 
random matrix theory can be used to predict the de- 
caying behavior of the fidelity M(t) in the golden rule 
regime, but not in the regime where the perturbation- 
independent decay of M(t) appears. The initial states, 
whose fidelity is found to have perturbation-independent 
exponential decay controlled by the Lyapunov exponent, 
such as narrow wave packets in the configuration space 
[3], coherent states [4], momentum eigenstatcs [9], etc., 
usually have wide spreading in the chaotic states \a) 
Such initial states, which are regular in the clas- 
sical limit, possess certain kind of coherence. It must be 
pointed out that the semiclassical theory does not pre- 
dict the perturbation-independent decay of the fidelity 
for an arbitrary initial state [14]. As shown in the previ- 
ous paragraph, the function Fm(E), equivalently, M(t), 
is basically determined by the quantity A a p in Eq. (9). 
In the golden rule regime, the LDOS Pl{E) is narrow, 
i.e., Cp a has many small components; as a result, for 
the quantity A a p, components of the initial states in the 
chaotic states, namely, (^o\(3){a\^o) , are not equally ef- 
fective in the whole (quasi)energy region of Hq and H. 
This should suppress the coherence of the initial states, 
which makes it possible to treat the components of the 
initial state taking part in the evolution of M(t) effec- 
tively as random numbers, like in the random matrix 
theory treatment in Refs. [4,6,14]. On the other hand, 
when the perturbation is so strong that the width of the 
LDOS is comparable to the (quasi)energy band width of 
Hq, the coherence possessed by the initial states should 
play a non-negligible role in the evolution of the fidelity. 
The random matrix theory can not be applied in this 
case, since it cannot describe the coherence in the initial 
states. 



III. NUMERICAL STUDY OF FIDELITY AND 
ITS FOURIER TRANSFORM 

A. The model 

The unperturbed Hamiltonian H of the kicked top 
model used in this paper is 

H ° = ^ s y + § s *Y, d ( t - nT )> ( 14 ) 

n 

where S is the total angular momentum, r is the period, 
and K is a parameter adjusting the strength of the kicks. 



Without loss of generality, the period r is set to be unit, 
r = 1 . The model describes a vector spin that undergoes 
a free precession around the y-axis and is periodically 
perturbed by kicks around the z-axis. The time evolu- 
tion of an initial state |$o) at t = is governed by the 
Floquet operator 

F =expH|^ 2 ]expH^], (15) 

where h has been set to be unit. The classical limit of the 
system, which is obtained by letting S — > oo with 1/S 
serving as the effective Planck constant, is fully chaotic 
for K > 9. 

The perturbation kV is chosen in the same way as in 
Rcf. [4], i.e., a slightly delayed periodic rotation of con- 
stant angle around the rr-axis, 

V=^-S x J2S(t-nr-e). (16) 

n 

The Floquet operator of the system H is 

F = cxp[-i^S x }F , (17) 
which gives the time evolution of an initial state, 

|*(t))=F t |$o>, (18) 

where t = 0,1,2,.... For this model, discussions in 
Sect. II are still valid, with the eigenstates of H and H 
changed to those of the Floquet operators F and F, re- 
spectively, with the corresponding eigenenergies changed 
to quasi-energies. 

In numerical calculations, for the sake of convenience, 
the perturbation parameter re is written as 

re = j x 1CT 3 . (19) 

Unless addressed explicitly, the value of S is 500 for the 
numerical calculations discussed in this paper. Since the 
time t in the evolution equation (18) takes integer values 
only, according to Eq. (11), the domain of the variable 
E in the function Fm{E), which is [— 47r,47r] in the gen- 
eral case, can be reduced to [— tt, tt], with the definition 
of Fm{E) changed accordingly, namely, 

F™ d (E)=J2 F M( E + 2n7r ^ ( 2 °) 

n 

where n takes the possible values in (0, ±1,±2). The 
function F$ d (E) is more closely related to properties of 
M(t) than the basic one in Eq. (12). It is this function 
that will be used in what follows in the numerical inves- 
tigation of the kicked top model and, for brevity, it will 
be denoted by F M (E). 

Initial states studied in this paper are in most cases 
coherent states of the SU(2) group [16,17], 
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\z) = Ae z " s+ \M = -S) (21) 

z = -tan-e**, (22) 
2 ' 

where \M) is the eigenvector of S z , with the eigenvalue 
M (M = -5, -5 + 1, . . . , 5), and A is the normaliza- 
tion coefficient. The coherence in a coherent state \z) 
can be easily seen when it is expanded in the states \M). 
However, in the study of the fidelity, \z) is required to 
be expanded in the eigenstates of the Floquet operators 
Fq and F. Although the expanding coefficients in such 
chaotic states must contain information on the coherence 
possessed by the coherent state, manifestation of the in- 
formation is not so easy. In Sect. HID, we will show nu- 
merically that the function Fm(E) introduced in Sect. II 
can supply some of the information. 



M 1 (t d ) = -J2M(t d ,$ i ), (23) 



we also calculate another averaged M(t), denoted by 
M2(td), by taking the logarithm of M(td) before per- 
forming the summation, 

1 - 

lnAfafa) = -2^1nM(t d) $i). (24) 
i=i 

Here M(t, i>i) is the fidelity of an initial coherent state 
defined in Eq. (1), with the dependence on writ- 
ten explicitly. 



Perturbation-independent decay of fidelity and 
quantum fluctuations 



Basic properties of the fidelity M(t) in the kicked top 
model have been studied in Ref. [4]. Here, as mentioned 
in the introduction, we are interested in a perturbation- 
independent decay of the averaged M(t) with initial co- 
herent states, which decays more slowly than the expo- 
nential decay predicted by the semiclassical theory. In 
order to understand this phenomenon, we study the indi- 
vidual M(t). In the regime of the perturbation parameter 
where the perturbation-independent decay appears, indi- 
vidual M{t) is found to have large fluctuations. See Fig. 1 
for some examples of K — 13.1 and j ' = 6, where M(t) 
are plotted in the logarithm scale. On the other hand, 
in the golden rule regime, e.g., K = 13.1 and j = 1, no 
obvious fluctuation in individual M{t) is observed. The 
large fluctuations are in fact quantum effects caused by 
the coherence in the coherent states. According to the 
arguments given in Sect. II, they should appear when 
the expanding coefficients of the states \(3) in |a) spread 
in the whole quasi-energy band of F$ effectively, that is, 
when the half width of the LDOS becomes comparable 
to the bandwidth of the quasi-energy, which is found nu- 
merically for j > 3. In the classical limit S — > oo, the 
fluctuations should disappear. 

As shown in Fig. 1, each M(t) has a 'shoulder' in a 
short initial time, denoted by t s in what follows, within 
which it changes slowly. The exponential-type decay of 
M{t) appears only after t s . In order to study the aver- 
age decaying behavior of M(t), their 'shoulders' should 
be subtracted. To this end, we shift the time variable 
t to td = t — t s and take average over M(i<j). Here we 
meet a problem, since individual M(t) has large fluctua- 
tions. The averaged M(td) may be mainly determined by 
a small fraction of the M(td) taken for averaging, if they 
have extraordinarily large values. Therefore, in addition 
to the standard averaged M(td), denoted by Mi(ta), 
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FIG. 1. The fidelity M(t) of three initial coherent states 
chosen arbitrarily, when K = 13.1 and j = 6. The dotted 
straight line shows the exponential decay with the correspond- 
ing Lyapunov exponent A = 1.65. 

Numerical results for M\(td) (upper) and M^td) (bot- 
tom) are shown in Fig. 2, with increasing k (j from 1 to 
8), when K — 13.1. In calculating the two quantities, 
t s for each M(t) of j = 1 and 1.5 was taken to be the 
first t at which M(t + 1) < 0.9. For j > 2, t s was de- 
termined by the first t at which M(t)/M(t + 1) > 1.5. 
Only states satisfying M(t s ) > 0.85 were used in aver- 
aging. The reason of taking two methods in calculating 
t s is that M(t) decays slowly at j = 1 and 1.5. The 
number TV of coherent states taken for averaging is 1000. 
Results in Fig. 2 for M\(td) are in consistence with those 
in Ref. [4], namely, the decay of M\(td) saturates at an 
exponential decay which is slower than the one with the 
Lyapunov exponent. However, when the large fluctua- 
tions in individual M(t) are suppressed by the second 
averaging procedure, Fig. 2 shows that the decay of the 
averaged fidelity M 2 (td) saturates at an exponential de- 
cay that is close to the Lyapunov decay. The deviation 
of the saturation decay of M\{td) from the one predicted 
by the semiclassical theory is caused by the large fluctu- 
ations in individual M(t). In fact, this can also be seen 
roughly in Fig. 1. 
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calculate their decaying exponents 71 and 72, 




1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 

2 4 6 8 10 12 14 



FIG. 2. Decay of averaged M(t) of initial coherent states, 
from the golden rule regime (j — 1) to the perturba- 
tion-independent regime, for K = 13.1. The difference 
between Mi and M2 lies in the averaging procedure (see 
Eqs. (23) and (24)). The dotted straight lines represent the 
exponential decay with the Lyapunov exponent A = 1.65. 



Mi oc exp(-7i* d ) M 2 (t d ) oc exp(-7 2 t ci ), (25) 

the values of which are presented in Fig. 3. Numeri- 
cally, 7i (i = 1, 2) are calculated by the best linear fitting 
to lnM^trf), with the fitting lines fixed to the values of 
lnMi(O) at td = 0. The points td used in fitting are 
those satisfying 1 > M 2 (td) > 0.008. In fig. 3, we see 
that, in the golden rule regime (j = 1 and 1.5), 71 is 
close to 72 because of the small fluctuation in individual 
M(t), and both of them are close to the half width of 
LDOS, r^, as predicted by the random matrix theory. 
The difference of the three quantities becomes obvious, 
when j > 2. In the parameter regime j > 6, 72 fluctu- 
ates around the Lyapunov exponent A = 16.5, while 71 
is obviously smaller than A. The value of 72 at j = 14 is 
obviously larger than the Lyapunov exponent. This may 
be due to the reason that the perturbation cannot be re- 
garded as small at this perturbation value. In fact, when 
the perturbation is so strong that the classical pertur- 
bation theory breaks down, M(t) may decay faster than 
the Lyapunov decay [6]. As discussed above, when S is 
increased and we are more and more close to the classical 
limit, the fluctuations in M(t) shown in Fig. 1 should be- 
come small and small, which would make 71 approach to 
7 2 . Indeed, we calculate the case of S — 1000 (k reduced 
to half to those of S — 500), where 72 are found to be 
fluctuating around A as well and 71 a little larger than 
the values at S = 500. However, 71 of S = 1000 are still 
obviously smaller than A, showing that S = 1000 is still 
not in the deep semiclassical regime. 
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FIG. 3. Values of 71 and 72 of the exponential decay of 
Mi(td) and M 2 (t d ), respectively, at different values of the per- 
turbation parameter, together with Tl, the half width the av- 
eraged LDOS, and Vf, the averaged half- width of the best fit- 
ting Lorentzian form to the central part of Pf(E) = |.Fm(£0|. 
The value of the Lyapunov exponent of the underlying classi- 
cal dynamics, A = 1.65, is indicated by the horizontal dotted 
line. 

To show the difference between the decay of Mi(td) 
and that of M2(td) in Fig. 2 in a quantitative way, we 



C. Dependence of fidelity on initial states 

In different regimes of the perturbation parameter, the 
initial state can play different roles in influencing the be- 
havior of the fidelity. In the Golden rule regime, the 
LDOS is narrow and, as discussed in Sect. II, the ran- 
dom matrix theory can be used to predict the behavior 
of the fidelity. In this regime, there should be no differ- 
ence between the decaying behavior of the fidelity of an 
initial coherent state and that of an arbitrary initial state, 
that is, the decay of the fidelity should be initial-state- 
independent. Indeed, in a study of the influence of sub- 
Planck scale structures [18] on the fidelity in Ref. [14], 
the averaged fidelity of the initial states 

\$ {T))=eM-iH T)\<f> c ), (26) 

where |$q) are coherent states, has been found numeri- 
cally to be independent on the parameter T and almost 
the same as that of initial coherent states. We compare 
the averaged fidelity of initial coherent states with that 
of initial random states of the form 

K andom > = E c «i Q >' ( 2? ) 
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where c a are random complex numbers satisfying the nor- 
malization condition. It is found that the former ones 
have 'shoulders', while the latter ones do not. With the 
'shoulders' of the former ones subtracted, the decay of 
the corresponding averaged fidelity is found to be close 
to each other in the golden rule regime. Some of the re- 
sults are plotted in Fig. 4, where 1000 initial states were 
used in averaging for each case. 
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FIG. 4. Comparison of the decay of the averaged fidelity 
of initial coherent states and that of initial random states in 
Eq. (27). 

In the perturbation-independent regime of the fidelity 
of initial coherent states, M(t) of |$o(T)) is found to be 
decaying in part of the evolution time in a way similar 
to that of initial coherent states [14]. This implies that, 
although the system H is classically chaotic, the time 
evolution operator exp(—iH n T) does not destroy the co- 
herence possessed by a coherent state. Indeed, we find 
that random superpositions of coherent states defined by 



i*o(«c))=x;. 



\ z n), 



(28) 



where (p n are random phases and \z n ) are coherent states 
in Eq. (21), still have some coherence properties of the 
coherent states. The averaged fidelity with $g s (n c )) be- 
ing initial states at the perturbation parameter j = 6 
are presented in Fig. 5. We see that the averaged M(t) 
with n c > 1 has a rapid decrease between t = and 
1, which is faster than the Lyapunov decay for n c > 8. 
This indicates that part of the coherence in the coher- 
ent components \z) of the initial states <I>Q S (n c )) is de- 
stroyed by the random phases in the initial states. In 
fact, when (a|$o)i the components of initial states in the 
eigenstates \a), can be treated as random numbers, the 
random matrix theory analysis predicts that the decay of 
the averaged M(t) is controlled by the averaged LDOS 
[14] , the half width of which is larger than the Lyapunov 
exponent at j = 6 (Fig. 3). 

After the initial transient rapid decrease, the rema- 
nent coherence in the initial states |$q (u c )) plays a role 



similar to that of initial coherent states. Indeed, the av- 
eraged M(t) changes slowly in some short time intervals 
after t = 1, which are analogues of the 'shoulders' of 
the fidelity of initial coherent states, then, decreases ex- 
ponentially in a way similar to that of initial coherent 
states, until they become close to the saturation value, 
which is about 1/(25*). The larger the value of n c is, the 
smaller the remanent coherence will be. The fidelity of 
random superpositions of \a) in Eq. (27) decays faster 
than the Lyapunov decay shown by the dotted straight 
line in Fig. 5, as predicted by the random matrix theory. 
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FIG. 5. The averaged fidelity with |$o s (n c )), a random su- 
perposition of coherent states, being the initial states. The 
average is taken in the same way as for M\(td)- The solid 
squares connected by the solid line show the averaged fidelity 
of initial random states in Eq. (27). The dotted straight line 
indicates the exponential decay with the Lyapunov exponent 
A. j = 6 for the perturbation parameter. 



D. Properties of Fm(E) and relation to properties of 

M(t) 

As shown in Sect. II, the Fourier transform of M(t), 
the function Fm(E), is directly determined by the states 
| a), \0) and |$o)- ln this subsection, we study properties 
of the function Fm(E) in the kicked top model numer- 
ically, in particular, its absolute value part pf(E) and 
its phase part 0p{E). Unless addressed explicitly, initial 
states |<E>o) studied in this subsection are coherent states. 
It will be shown that properties of the central part of 
Fm{E) are closely related to the decaying rate and the 
width of the 'shoulder' of the averaged fidelity M(t). 

Typical shapes of pf(E) and 9p(E) in cases of dif- 
ferent values of the perturbation parameter are shown 
in Fig. 6, where the same coherent state is used. For 
Pf(E) of j = 1 and 2, the values of pf are plotted in the 
logarithm scale. Each pp(E), symmetric with respect 
to E — 0, is composed of three or more parts, with a 
main central part. A manifestation of the coherence in 
the coherent state |$o) is that the corresponding phase 
functions 8f(E) have related divisions, as can be seen in 
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Fig. 6. (If |$ ) is taken as a random state |$r andom ) m 
Eq. (27), values of the function 6p(E) will be close zero.) 
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FIG. 6. Typical shapes of Pf(E) and 0f(E) of the same 
coherent state, for j — 1, 2, 3, 4. Pf{E) are plotted in the 
logarithm scale for j — 1 and 2. 
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FIG. 7. Lorentzian fit to the central part of the Pf(E) in 
the j = 4 case in Fig. 6. 

Numerical results for the phase function Of{E) show 
that it is approximately linear in the region of E corre- 
sponding to the central part of pp (see Fig. 6 for exam- 
ples), 

e F {E) « kgE, (29) 

where kg is the slope of the approximate linear behavior. 
Distributions of ke are shown in Fig. 8 for some pertur- 
bation parameters. The positions of the peaks of the 
distributions are close to kg = 1, while the distribution 
of j = 4 has another small peak at kg — 1.5. To illustrate 
the influence of 9f(E) on behaviors of M(t) clearly, let 
us write M(t) in the form 

M(t) = f p F {E)e Wt{E) dE 1 (30) 



The central part of each pp is found to be fitted well 
by the Lorentzian form, with a width denoted by Tp, to 
be distinguished from Yl for the width of the LDOS. See 
Fig. 7 for an example of the Lorentzian fit to the central 
part of the pp in Fig. 6 of j=4. The Lorentzian form used 
in fitting is pl{E) + a, where Pl{E) is the function on 
the right hand side of Eq. (3) and a is a fitting constant. 
Therefore, the closeness of the central part of pf(E) to 
the fitting Lorentzian form does not predict directly an 
approximate exponential decay of M(t) with an expo- 
nent Tj?. In fact, the contribution of the central part of 
Fm{E) to M{t) in Eq. (11) gives both positive and nega- 
tive results. The final positive results for M (t) come out, 
only when the contributions of the other parts of Fm are 
taken into account. However, the average values of Tp, 
which are shown by solid squares in Fig. 3, are found 
numerically to be close to the values of 71 for Mi(td), in 
the whole parameter regime of j from 1 to 14. When 5* 
is taken equal to 1000, the change of the corresponding 
results of Tp are found to be smaller than that of 71. 
Since 71 increases only a little, when S is changed from 
500 to 1000, it is not clear at the present stage whether 
or not can approach to 71 in the classical limit. 



where 8t(E) = 9p(E) — Et. Some examples of the be- 
havior of t (E) are shown in Fig. 9. With increasing t, 
the approximate slope of 9t{E) in the central region of 
Fm{E) changes from positive to negative. For t > 2, the 
larger t is, the steeper 9t(E) will be. As shown below, 
this feature of t (E) can be used to estimate the width 
of the 'shoulder' of M(t). 
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FIG. 8. Distributions of kg, the slope of the approximate 
linear behavior of 8f{E) in the central region of Fm(E). 
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FIG. 9. Change of 8 t (E) with t, for the F (E) of j = 4 in 
Fig. 6. 
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In order to give an estimation of the 'shoulder' width of 
M(t), we note that the exponential decay of M(t) should 
appear when 9 t {E) is steep enough, specifically, when 
6t(— Wf/2) is comparable to ir, where Wf is the width 
of the region of E occupied by the main central part of 
Pf- Quantitatively, we use t s g to denote the average of 
the largest t satisfying 



Wf 

\ke-t\— < R c <ir, 



(31) 



where R c is a quantity used to show the closeness to n. 
The average width of the 'shoulders' of M(t), denoted 
by t S M, are calculated from the averaged M(t) directly. 
Since t has integer values only in the fidelity M(t) of the 
kicked top model, we use the linear interpolation in cal- 
culating t S M by M(t S M) = M c , where M c is a quantity 
measuring the lower border of the 'shoulder' of the aver- 
aged M(t). In studying the relation between t S M and t s g, 
for a given value of M c , the value of R c is determined by 
the smallest value of 



(t s g ~ t s MY 



(32) 



in the parameter regime 1 < j < 8. Figure 10 shows two 
examples of M c = 0.9 with R c ~ 0.75 and M c = 0.8 with 
R c ~ 0.9, where we see that t s g w t S M for j between 1 and 
8. The difference between t s g and t S M becomes increas- 
ingly large, when j is larger than 8, which is related to 
the fact that the perturbation kV is not quite small for 
j > 8. Variation of R c with respect to M c is presented 
in Fig. 11, where R c is seen to be almost linear with M c 
for M r between 0.7 and 0.95. 



FIG. 10. Two examples of the relationship between t S M 
and t s 0. They are close to each other, when j is not too large, 
namely, 1 < j < 8. 
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FIG. 11. Dependence of R c on M c , required by the smallest 
value of (32). 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we study the crossover of the fi- 
delity from the golden rule regime to the perturbation- 
independent decay regime. We show that the deviation 
of the perturbation-independent decay of the averaged 
fidelity from the prediction of the semiclassical theory 
is due to large quantum fluctuations in individual fi- 
delity caused by the coherence in the initial coherent 
states. When the quantum fluctuations are suppressed 
by an appropriate averaging procedure, the perturbation- 
independent decay of the fidelity is found to be close to 
the semiclassical prediction. 

We also show that when the effective Planck con- 
stant is not small enough to guarantee the validity of 
the semiclassical theory, the crossover does not occur at 
a point of the perturbation parameter, but takes place 
in a region of the parameter (Fig. 3). The crossover 
in this case can be analyzed by neither the semiclassi- 
cal theory, because the effective Planck constant is not 
small enough, nor the random matrix theory, since the 
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perturbation-independent decay of the fidelity is initial- 
state-dependent . 

Furthermore, a function Fm(E) that is in fact the 
Fourier transform of the fidelity, is introduced to study 
the fidelity in the crossover region. This function is de- 
termined directly by the initial state and the cigenstates 
of the two Floquet operators concerned. Numerically, 
for initial coherent states, the absolute value part and 
the phase part of the function Fm(E) are found to have 
related division into several sub-parts, which is a manifes- 
tation of the coherence possessed by the initial coherent 
states. In the whole crossover region, properties of the 
central part of Fm{E) are found to be closely related to 
both the exponent of the exponential decay of the av- 
eraged fidelity and the length of the short initial time 
period, within which the fidelity almost does not change 
due to the coherence possessed by the initial coherent 
state and the smallness of the perturbation in the classi- 
cal limit. 

The function Fm{E) is studied mainly numerically in 
this paper. A detailed analytical analysis of the function, 
e.g., by making use of its expression as the correlation 
function of f m (e) given in Sect. II, would supply further 
understanding of the behaviors of the fidelity, in analogue 
to the relation between the LDOS and the survival prob- 
ability. Information of the coherence in a coherent state 
may be hidden in its expanding coefficients in the eigen- 
states of classically chaotic systems. Numerical results in 
this paper show that the function Fm supplies a useful 
method of extracting such information. We would like 
to mention that some other quantities discussed in this 
paper may serve this purpose as well in other situations, 
e.g., the quantity A a p introduced in Sect. II, which is ir- 
relevant to the arbitrary phases that can be given to the 
initial coherent state and the quasi-energy eigenstates of 
the two classically chaotic systems. 
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